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Abstract 

We consider a type of evolution on {0, 1}™ which occurs in discrete steps whereby 
at each step, we replace every occurrence of the substring "01" by "10" . After at most 
n — 1 steps we will reach a string of the form 11 • • • 1100 . . . 00, which we will call a 
"stabilized" string. If we choose each bit of the string independently to be a 1 with 
probability p and a with probability 1 — p, then the time of stabilization of a string in 
{0, 1}™ is a random variable with values in {0, 1, ... n — 1}. We study the asymptotic 
behavior of this random variable as n — > oo and we determine its limit distribution after 
suitable centering by mean and scaling by variance. The case P = \ turns out to be 
different from any other values of p; in this case the limit distribution is %3, which is 
distribution of the length of (Z\, Z%, Z3) g R 3 where Z\, Z2, Z3 are independent -/V(0, 1) 
random variables. For p 7^ |, the limit distribution is Gaussian N(0,p(l — p)). The 
tools used in our arguments are a natural interpretation of strings in {0, 1}" as Young 
diagrams, and a connection with the known distribution for the maximal height of a 
Brownian path on [0, 1]. 

1 Introduction 

For n G N and p G (0, 1) let denote the probability space consisting of strings in 
{0, l} n , where elements uj are assigned the product probability induced by assigning each 
bit independently to be a 1 with probability p or a with probability 1 — p. We consider 
the following kind of 'evolution' for strings in f2^: for a fixed cu G 0^ we look at all the 
occurrences of the substring "01" in oj and we replace every one of them by "10". By doing 
so, new instances of "01" may be created (for instance 0101 1— > 1010 creates a "01" in the 
middle). We repeat this process until we reach a string of the form 11 • • • 1100 • • • 00. 

A concrete example: say we have n = 8 and we start with the string u = 01101011. 
Then our evolution produces the following strings before stabilizing: 



01101011 H- 10110101 ^ 11011010 ^ 11101100 i-> 11110100 i-> 11111000, 

and stabilizes after 5 iterations because there are no more instances of "01" to be found. 

This evolution has a cute interpretation, the original inspiration for the problem, as a 
line of confused soldiers, see [5]. There is also an interpretation as particles whose motion 
is restricted to one dimension. Imagine that each 1 in the string is a particle that would 
like to move to the left-hand side of the string and that each is an empty space. At every 
iteration, if a particle has an empty space to its left, it will move into that space. This is 
exactly the replacement rule 01 — > 10. Equally well, one can think of the 0's as the particles 
which are trying to move as far right as they can and the l's as open spaces. The process 
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will stabilize when all of the particles have moved as far as they can go. In this guise, 
the problem is a kind of deterministic analogue of certain "exclusion processes" (see for 
instance [3] Chapter 8) except in our case the initial condition is random, but the evolution 
is deterministic. 

It is elementary to show that this process must stabilize after at most n — 1 iterations, 
and that afterwards we will obtain a string of the form 11 • • • 1100 • • • 00. The number of 
iterations until we reach such a final configuration is a random variable on the probability 
space Qn . We will call this random variable the stabilization time and denote it by T% : — > 
{0, 1, . . . , n — 1}. Going back to our concrete example above, we have Tg (01101011) = 5. 

In this paper we will examine the limit distribution for the random variable T% in the 
limit n — > oo and for varying values of p. One of the points of interest is the fact that the 
limit distribution in the case p = ^ is qualitatively different than it is for other values of p. 

Theorem 1.1. We have the following weak limits for the distribution of the random variable 
TS in the limit n — > oo: 



In the case p < ^: 
In the case p = \ : 



Tn 2 ^ 



In the case p > 



I . 

2 ■ 



T'n -pn 



X3 



n 



N(0,p(l-p)) 



In the above limits, N (0,p(l — p)) is a mean zero Gaussian variable with variance p(l —p), 
and X3 ~ \J Z\ + -Zf + is the Euclidean norm of a vector of three independent standard 
N(0, 1) Gaussian variables. This has density: 

/ \ 8\/2 n — 2t 2 r r r\ 

dX3\ x ) = — i= x e dx for x > 
vvr 



The proof of the theorem comes through connecting several ideas from combinatorics 
and probability theory. First, we notice a natural connection between strings in On and 
Young diagrams. Using Young diagrams as a tool, we can analyze the special case of strings 
that begin with a and end with a 1 and relate the stabilization time to a simple random 
walk in one dimension. Using this connection to the random walk, we can find the limit 
distribution for the stabilization time in the special case mentioned above using the central 
limit theorem and Donsker's theorem. Finally, we show that the limit distribution for the 
special case of strings is actually the same as the limit distribution for general strings. 

We have divided the proof into three lemmas which are stated below and each discussed 
and proved in their own sections. 

Lemma 1.2. Given a string uj G {0, l} n , there is a a natural ID random walk associated 
with uj which takes a step up for every in uj and takes a step down for every 1 in uj. To be 
precise, for < k < n let Sf. = 2? = i (1 ~ 2wj) . In the special case where uj\ = and ui n = 1, 
we have the following explicit relationship between the random walk Sk and the stabilization 
time Tn{uj): 
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T n » = ^+ max S k -^-l 

Z l<k<n Z 

Remark 1.3. The proof of this lemma comes by mapping each string to a Young diagram 
in a natural way. The stabilization time of string turns out to be equal to the property 
called the depth of the Young diagram. In the special case uj\ = and uj n = 1, the depth of 
the Young diagram is also found to be equal to the above expression for our random walk. 
The proof of this lemma is discussed in Section [2j 

Lemma 1.4. Let Qn = {u> 6 {0, l} n : uj\ = 0, uj n = 1} be the probability space where each 
bit, except for the first and last, is chosen independently at random to be a 1 with probability 
p and to be with probability 1 — p. The stabilization time of a string is a random variable 
on this probability space which we will denote by T% ■ f^n — )■ {1, 2, . . . , n — 1}. T% has the 
following weak limits as n — >• oo: 



In the case p < 



In the case p 



i 



1 . 

2 ' 



Tn - n(l - p) 



n 



JV(0,p(l-p)) 



In the case p > \. 



'"' ->JV(0,p(l-p)) 



'n 

Remark 1.5. Once Tn{u) = § + max i<yt<n — — 1 is established in Lemma 11.21 the 
proof of Lemma 11.41 is an exercise using the central limit theorem and Donsker's theorem. 
For p < \ some simple analysis shows that the term maxi<fc< n Sk does not contribute to 
the limit distribution, and the central limit theorem is enough to prove the limit. The case 
p > \ follows by symmetry between the 0's and l's in the problem. In the case of p = ^, 
the random variable X3 arises from a connection to Brownian motion. More precisely, if Bt 
is a Brownian motion and Mt = max s <j B s is its running maximum then: 

Mi - l -B x I X 3 

and for a simple random walk Sk, we have by Donsker's theorem that: 

inaxi<Kn S& — hSm 1 
In 2 

The full details of the proof of this lemma are displayed in Section 3. 

Lemma 1.6. Let fir be the law of the special case random variable Tfi and let $i be the law 
of the random variable Tn. Let b~§ be the unit mass at 0. 
For p = \ these measures are related by: 

| n + 1 r-^rt — r— 1 _| 

r=0 
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For p 7^ \ these measures are related by: 

^ (l_ p ) p n-r_ p(1 _ p) n-r 

^ = 2^1 5 ° + ^ 2^1 

^ r=0 ^ 

Moreover, from this relation, we can show that the random variable T% will have the 
same limit distribution as the special case random variable T^. 

Remark 1.7. The formulas stated in Lemma 1.6 are obtained by carefully counting the 
number of ways of adding leading l's on the left and trailing O's on the right of a string of 
special type Cj £ fin, and by exploiting the obvious fact that such additions of leading l's 
and trailing O's does not affect stabilization time. Hence, the stabilization time for a general 
string to is equal to the stabilization time for the special case substring Cj one obtains from 
u> by deleting any leading l's or trailing O's. This gives the relationship between \x and fi in 
the lemma. One can see that this is a convex combination of the probability measures jl r , 
which is heavily weighted towards higher values of r. Some elementary analysis is used to 
show that this weighting is such that [i and jl have the same limit distribution. The details 
of this lemma are given in Section 4. 



2 Connection to Random Walk in the Special Case 

In this section we aim to prove Lemma 11.21 which gives the following explicit formula for 
the stabilization time of w G {0, l} n satisfying u\ = and u n = 1. 

T n » = ^+ max S k -^-l 

Z l<k<n Z 

This formula, once the statement is known, can be proved directly by induction. How- 
ever, the proof by induction obscures the mechanics of what is really happening in the 
evolution. A more illuminating solution, which is how the formula was first discovered, is to 
map the strings into Young diagrams in a particular way. When viewed as Young diagrams 
in this way, the evolution process becomes particularly simple, and then the explicit formula 
can be seen more directly. This is the proof we will present in this section. We will begin 
by creating a map from {0, l} n to the space of Young diagrams. 

Definition 2.1. Let Y denote the collection of all finite sets Y C N x N which have the 
following property. 

(i,j) £ Y=J> (i',f) e Y for all i',f £ N such that i! < i and f < j. 

The empty set is also counted in Y. It is customary to represent a set Y € Y as a picture 
containing a "collection of boxes" (square boxes of side 1), where for every £ Y we 
take in our picture the box which has (i, j) as its top right corner. In this guise, the set Y is 
called a Young diagram, see e.g. [4J. The following example illustrates this representation. 



Y = {(1, 1), (2, 1), (3, 1), (4, 1), (1, 2), (2, 2), (1, 3), (2, 3)} = 
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Definition 2.2. For every uj G {0, 1}™ with uj\ = and u; n = 1, let [7 = be the 

number of l's in the string uj. Let B be a J7 x (n — U) grid on the x-y plane. Reading a; 
from left to right, we construct a path starting from the top- left corner of B by drawing a 
line horizontally to the right whenever we encounter a in a; and a line vertically downward 
whenever we encounter a 1 in uj. Since there are U "l"'s and n — U "0"'s to be found in oj 
we will get a path from the top-left corner in B to the bottom-right corner in B. Here is 
an example of the path generated by the string uj = (0, 1, 1, 0, 1, 0, 1, 1): 



Define tt(uj) C N x N to be the unique path constructed in this manner. In our example, 
tt(0, 1, 1, 0, 1,0, 1, 1) = {(0, 5), (1, 5), (1, 4), (1, 3), (2, 3), (2, 2), (3, 2), (3, 1), (3, 0)}. The set of 
boxes under this path defines a Young diagram, which we denote Y{uS). In our example, 
this is the following diagram. 



Y(u) 



Definition 2.3. Let Y G Y be a Young diagram and let G Y. We say that is an 
exposed corner of Y if (i + G" Y and + 1) ^LY . The corner cutting map K : Y — > Y 
is the map that removes all the exposed corners of a Young diagram as follows. 

K(Y) :=Y\ G Y \ is an exposed corner of Y} 

It is easily verified by the definition that removing the exposed corners does indeed yield 
another Young diagram. That is to say, K is a well-defined function from Y — > Y. 

Proposition 2.4. Let uj G {0, 1}™ be a string and let uj' be the string obtained from uj after 
one stage of the evolution, that is after replacing instances of "01" with "10" once. Then 
the Young diagram Y(uj') is obtained by cutting the corners of the Young diagram Y(u>). 

Y{u')=K{Y{uj)) 



Proof. Consider the path tt(uj) defined in Definition 12.21 Any instances of "01" in ui will 
correspond to a horizontal segment followed by a vertical segment, while instances of "10" 
correspond to a vertical segment followed by a horizontal segment. As such, the evolution 
"01" —> "10 will translate into the following pictorial evolution for the path it(uj). 



These instances of "01" are correspond exactly to the exposed corners of Y{oj) since they 
have no neighbors above them or to their right. We can see that this evolution is precisely 
removing these exposed corners of Y(oo). ■ 
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Definition 2.5. Let Y G Y, Y 7^ be a non-empty Young diagram. We define: 



Depth(Y) = maxji +j- G Y} 

We also set the convention Depth(0) = 0. 

Proposition 2.6. Let Y G Y, Y ^ 0. T/ien Depth(K(Y)) = Depth(Y) - 1. 

Proof. Suppose G Y such that i + j — 1 = Depth(Y). Clearly (i + l,j) G" Y and 

(i, j + 1) Y, since assuming otherwise contradicts i + j — 1 = Depth(Y) is maximal. Then 
is an exposed corner, so (i, j) G" K(Y). Hence Depth (If (Y)) < Depth(Y). 
Now, if either (i - l,j) G K(Y) or (i,j - 1) G K(Y) then Depth(if(Y)) >i + j -2 = 
Depth(Y) — 1. By combining the inequalities, we conclude Depth(if (Y)) = Depth(Y) — 1. 
Otherwise, both (i — 1, j) ^ K(Y) and — 1) ^ if (Y). But this only happens in the case 
Y = {(1, 1)}, in which case K(Y) = and the proposition holds by Depth(0) =0. ■ 

Lemma 2.7. Let uj G {0, l} n and let T(uj) be its stabilization time. Then we have the 
following relation. 

T(uj) = Depth(Y(uj)) 

Proof. The proof follows by induction on Depth(Y(w) using the previous two propositions. 
For the base case, if Depth(Y(w)) = 0, then Y = so uj is already stable and the stabilization 
time is and the result holds. Now assume that for all uj with Depth(Y(cj)) = k — 1, 
T(uj) = k — 1. Given any uj with Depth(Y(w)) = k, let uj' be the evolution of uj by one 
step. Then, by the last two propositions we have that Depth(Y(u/)) = Depth(if (Y(lj))) = 
Depth(Y(cj)) — 1 = k— 1, and so by the induction hypothesis, T(uj') = k — 1. Since we have 
applied one step to get from uj to a/, T(uj) = T{uj') + 1 = k. ■ 

Proposition 2.8. Let uj G {0, l} n , and for < k < n let S k = Yh=i0- ~ 2uj i)- Let 
U = U(uj) be the number of 1 's in the string uj. In the case that oj% = and uj n = 1 we have 
the following explicit relationship between Sk and the depth of the Young diagram Y(uj): 

Depth(Y(uj)) = U+ max S k - 1 

0<fc<n 

Proof. Let tt(uj) be the path associated to uj in Definition 12.21 An elementary computation 
shows that S is related to tt(uj) by the following formula: 

{(k,S k ):0<k<n} = {(i + (U-j),i-(U-j)) : (i,j) G tt(uj)}. 

This formula merely says that (because of the specifics of how each of S and vr(w) are 
constructed from uj) the set of lattice points S(uj) := {(k,S k ) ■ < k < n} is obtained out 
of it(uj) via 45° degree rotation and dilation by The verification of the formula is left 
as exercise to the reader. An illustration of how tt(oj) and S(uj) look in a concrete case is 
shown in the next picture, drawn for uj = (0,0, 1, 1,0, 1, 1, 1,0, 1) G {0, l} 10 . 



Now, if uj\ = and uj n = 1, we know that the boundary of the Young diagram Y(oj) is 
precisely the path tt(uj). (The fact that ui\ = and uj n = 1 is needed here because otherwise 
there are some points in tt(uj) which are not included in Y(lo)). Since the Depth(Y(u;)) is 
achieved somewhere on its boundary, we have that: 

Depth(Y(w)) = maxjj + j - 1 : G Y(u)} 

= max{i + j - 1 : (i, j) G vr(o;)} 

= U + max {i — (U — j) : (i, j) G 7r(w)} — 1 

= [/ + max{5 fc : < fc < n} - 1 

The last equality follows from the map between S and tv(lj) described above, and the 
complete equality is our desired result. 

■ 

From here the result of Lemma 11.21 is a simple corollary obtained by using the elementary 
fact that U = \ (n — S n ). 



3 Limit Distribution in the Special Case 

In this section we aim to prove the aforementioned weak limits for the random variable 
Tn, the stabilization time for strings from the probability space Qn = G {0, 1}™ : U)\ = 
0, u) n = 1}. We will use the result from Lemma 11.21 connecting T% to the random walk 
associated with u, Sk = ^2i=i0- ~ 2wi)- Recall that: 

r£(u,) = ^+ max S k -^-l 

Z l<k<n Z 

Lemma 3.1. Let Xi,X2, ■ ■ ■ ,X n _2 be i.i.d random variables which take the value —1 with 
probability p and 1 with probability 1 — p. For < k < n — 2, let W k = ^f=i be the 
random walk which takes the X, L 's as its steps. Then: 

T£ = - + max W k — 

2 0<fc<n-2 2 

Proof. Make the identification that Xi = (1 — 2cjj+i). Since u\ = 0, co n = 1, we have that 

5fc + i = 1 + Wk for < k < n — 2 and S n = 1 + W n -2 — 1. From this, it is clear that 
maxi</ c < n 

S k = l + maxo<fc< n _2 VFfc, and the result is then immediate from Lemma IL2J B 

Since Wk is the sum of many i.i.d. random variables, we are in a position to use tools 
like the central limit theorem and Donsker's theorem to find the limit distribution. We 
divide the remaining results into the case when p and the case when p =\. 

Proposition 3.2. For p > ^, we have that, as n — >■ oo: 

max < fc <n- 2 W k 



ii 



Proof. maxo<fc<oo Wk is the maximum height achieved at any time by a weighted random 
walk, which takes steps upwards with probability 1 — p < i and steps downwards with 
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probability p > |. It is a result from elementary probability that this is distributed like a 
geometric random variable related to the parameter q = 1 ^ l 2 p P ^ ^ < L 

max Wk = Geom(l — q) — 1 

Here, q is the probability that if the random walk Wk reaches height x, the walk will ever 
reach the height x + 1. The random walk can either reach x + 1 immediately by taking a 
step upward at this time step (with probability 1 — p) or it can step down now, and then 
take two steps upward eventually (these are independent, so this happens with probability 
p ■ q • q). Hence q satisfies the quadratic equation q = (1 — p) + pq 2 , the solving of which 
admits the stated value for q. To see that the maximum is indeed distributed as a geometric 
random variable, write: 



P(maxW fc = m) = P(3k : W k = l\3i : W = 0) • P(3k : W k = 2\3i : W* = 1) • . . . 
.... p(Vfc, W k ^ m + l\3i : Wi = m) 
= q-q...q-(l-q) = q m (l-q) 

from which the relationship to the geometric random variable is clear. In particular, 
E(maxH / fc) = j— ^ — 1 < oo and so we have by the Markov inequality that: 

/ maxQ</j< n _2 Wk \ < E (max < fc <oo W k ) 



as n — > oo 



Lemma 3.3. Let Cl^ = {uj G {0, l} n : Wi = 0,o; n = 1} 5e a probability space and Tn be the 
stabilization time on this probability space both as defined in Lemma \l.4\ For p > \,Tn has 
the following weak convergence to a Gaussian as n — )• oo; 

V "~"" ^iV(0,p(l-p)) 



in 

Proof. We have, using the results of Lemma 13. II and Proposition 13. 2\ that: 



np 



f (2p-l) 



max <K n _ 2 H/ fc - %^ ) - f (2p - 1 



max <fc<„_ 2 W fc 1 Wn-2 - "-(1 - 2p) 



+ ^p(i-p) 



Since max °-^~ 2 — =^ by Propositon 13.21 and because W n _2 is the sum of the i.i.d. 
random variables Xj with mean 1 — 2p and variance 4p(l — p), so by the Central Limit 
Theorem, we have weak convergence to a Gaussian: ^sr^M N(0, 4p(l - p)). ■ 
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Corollary 3.4. The same limit distribution as in Lemma 3.3 holds for p < ^. 

Proof. As described in the introduction, we can imagine the evolution process from two 
points of view: we can think of each 1 in the string is a particle trying to move to the left 
of the string, or we could equally well imagine that each is a particle trying to move to 
the right of the string. Each of these ways of thinking about the problem is equally valid, 
and gives the correspondence between p and 1 — p. 

More formally, one can identify the probability space Qn with the probability space &n~ p 
with the following bijection: interchanging the O's and l's and then reverse the order of the 
string. This bijection preserves the evolution operation of "01" — > "10" and hence shows 
that that Tn and Tn~ p have the same distribution for every n. ■ 

We now state "Donsker's Theorem" , which is used to find the limit distribution in the 
case p =\. - see e.g. Section 8 of [lj. 

Theorem 3.5. (Donsker's Theorem.) Let X\, X2, X3, . . . be i.i.d. random variables that 
are ±1 each with probability |. Let W n = ^2i<j< n Xj be the sum up to time n. (This 
is a simple symmetric ID random walk). For each n, define a piecewise linear function 
L n (t) that is linear on each interval that L n (—) = ^i. Then, we have weak 

convergence of the law of L n {-) to the law of a Brownian motion Bt,t £ [0, 1] in the sense 
of weak convergence in the space of continuous functions C[0, 1] equipped with the sup-norm 
metric: 

L»(.) => B { . } 

Lemma 3.6. Let Bt,t £ [0, 1] be a Brownian motion and let Mi = max s <( B s be it's running 
maximum. Then, in the case that p = \ we have the following weak limit: 

1 

rp2 _ n 1 

^^^m 1 -\b 1 

n I 



Proof. We know that T^ = ^ + max <fc< n _2 — Wn 2 2 , so it suffices to show that: 

m&x < k < n W k - \ W n 1 

= Mi - -Bi 

n I 



This is a consequence of Donsker's theorem. Indeed, let h : C[0, 1] — > K by h(f(-)) = 
su Pte[o,i] /(*) — 5/(1)- This is a continuous function on C[0, 1] with the sup norm (indeed, 
both sup te [ 0] i] f(t) = 11/11 and evaluation at time 1 are continuous in this space). Now, since 
h is continuous, we know that it respects weak convergence. Applying our h to the limit 
in Donsker's theorem, L n (-) =^> B^ gives h(L n (-)) =^ h(B^). By our definition of h and of 
L n , this is precisely that: 

h(L n (-)) = sup L»(t)-h»(l) 
te[o,i] 1 

max <fc<„ W k -\W n 



jn 
KB { .)) 

sup B t -]-Bi 

te[o,i] 1 

Mi - -B x 
2 
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Lemma 3.7. M\ — \B\ = X3 with probability density function: 

dx3 = -j= x 2 exp (— 2x 2 ) dx 



7T 

Proof. We verify this by computing the density of M\ — \B\. This is just a computation 
using the joint density for Brownian motion and its Maximum, which is readily calculated 
using the reflection principle - see for example [2] pg 95. The joint density function for 
Brownian motion and its maximum is: 

n 2 (26 -a) ( {2b-a) 2 \ 
p(M T = b,B T = a)= K 3 exp - v ' 

for b > a, b > 0, and otherwise. Now to get the density function for M\ — \B\, one just 
integrates the joint density for B t and M t along a line: 



2;r 



p(H = x) = j p (Mi = | + x, B x = y) dy 



-2x 

- [ i(y + 2x ) - y) ex P ( - — + it - — ~ I d y 



7T J ' r \ 2 

-2x V 

-2x 

8 ^ 2 exp (-2x 2 ) 



7T 



Combining all the results from this section gives the result of Lemma II. 41 

4 Limit Distribution in the General Case 

So far, we have results in the special case that the string u has u)\ = and uj n = 1. In 
this section, we will bootstrap off of these results to see that we have the same limit in the 
general case where there is no restriction on u)\ or cj n . 

Lemma 4.1. Let be the law of the special case random variable Tfi and let p^ be the law 
of the random variable T%. For p ^ i the measures are related by: 

p n+l _ (1 _ p) n+l ^ {1 _ p)p n-r _ p{1 _ p) n-r 

^ = 2p-^l 60 + ^ 2p-Tl 
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Proof. We begin by splitting the space {0, 1}" into disjoint subsets. For stable uj £ {0, l} n , 
all "l"s in uj must lie to the left of all "0"s in uj. As such, there are n + 1 stable strings in 
{0, l} n for which the time to stabilization is zero. For each 1 < i < n, there is precisely one 
such string with exactly i "l"s. These strings contribute the following value to fin. 

YsAl-PT^o = P ~ 2 l_ l P) So (4.1) 

i=0 p 
For non-stable uj E {0, l} n , we will introduce r, the number of elements of uj lying between 
the first "0" and the last "1", so that we may write uj = (1, 1, 0, xi, x r , 1, 0, 0), 
with x G {0, l} r . The time until stabilization of w is the time until stabilization of 
(0, xi, X2—x r , 1). For given r, the distribution of these times is precisely distributed like 
Tf + 2 because this string is in the special case. It remains to count the number of ways we 
can get r to be a fixed value. Outside of (x\,X2---x r ), we require one "1", one "0", and 
n — 2 — r instances of either "1" or "0" in arrangements that are unique for a given number 
of "l"s and "0"s. These strings contribute the following value to fi n for a given r: 

n x> m (i -PY— i ti + 2 = {1 ~ p) ^~'~_f~ p)n ~ r K + 2 (4.2) 

i=0 p 
By summing over all possible values of r, we get the desired result. ■ 

Lemma 4.2. When p = \, we have: 

i_n + l y-^n — r — 1 _i 

r=0 

Proof. The proof is the same as above. The only difference is that in the case p = i the 
sums we need to evaluate are arithmetic sums, instead of geometric ones. One can also 
obtain the result by taking the limit p — > \ and computing the indeterminate limits. ■ 

Notation 4.3. We introduce the notation Cr,n, a set of positive coefficients defined by the 
results of Lemma 14.11 and 14.21 so that the following holds true for all values of p: 

n 

r=2 

The reason for doing this is to unify the notation from the cases P = \ and p ^ i. To 
prove the result we are after, we only need three properties of the coefficients Cr, n - These 
are proven in the next lemma. 

Lemma 4.4. For each n, we choose a parameter 6 n , for which 9 n — > oo as n — > oo and 
^= — > as n — > oo. We will use 6 n = \n a \ where < a < | as an illustration of these 
properties. Given 8 n as defined above, we have the following properties of Cf n : 



1- ^O.n + ^r,n 
r=2 
n—6 n 



r=2 

n 



Env 

r=n—9 n +l 
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Proof. Property 1 holds since the Cr,nS arose in Lemma 14,11 and Lemma 14.21 as the division 
of the probability space into disjoint pieces. This can be verified directly by summing 
over the values for Cf jn . Once 1 is established, 2 and 3 are equivalent. To establish 2, we 
use the definitions of Cr, n which are found in Lemma 14.11 and Lemma 14.21 For p < | , we 
have the following equation. 



nv /x^nv P n+1 ~ ( 1 ~ P) n+1 , V^" (1 ~ P)P n ~ r+2 ~ PQ- ~ P) n ~ T+2 

r=2 ^ r=2 v 

~ 2p-l 2p - 1 ^ p 

r.n+1 I r.On+2 n n+l 

P \ P P n— »oo 

2p- 1 

The last line holds since 9 n — > oo as n — > oo. The proof is similar when p > \. Finally, 
when p = i, we have the following equation. 

n—On . 1 n—0 n . 

n + 1 n — r + 1 



r=2 r=2 



n + l 1 9 n + 2 n + l\ „_>oo 



Proposition 4.5. T/ie random variable Tn has the same limit distribution as that ofT%. 
That is: 

I*-n(pV(l-p))^ 



n 



where p is the limit distribution of Tn P ^ as established in Lemma \l-4\ 

Proof. For the whole proof, we fix a value of p. We will discard the superscript p in Cr,m P>n 
and p r in the remainder of the proof to simplify the notation. We assume that p > \ so 
that the limit in question is Tn ^ n . In the case p < \ the limit in question is Tn ~^ p ^ n 
instead, but the argument is identical. 

We use the following characterization of weak convergence for distributions p n : p n =>• p 
if and only if p n ((— oo, a}) — > /?((— oo, a]) for every a G M with p({a}) = 0. This is true 
since the semi-infinite intervals {(— oo,a] : a G R} are a so-called convergence- determining 
class. See pQ example 2.3 pg 18. 

Let A = (— oo, a], a G M be an arbitrary semi-infinite interval. It will be useful in our 
study here to use the notation yJnA + pn = ^fn (— oo, a]-\-pn = (— oo, ^fna + pn] . Our aim 
is to show that: 



P 



(^i-J^eA\-p{A) = \p(T n eV^A + pn)-p(A)\ 

= \(p n ) (y/nA +pn) - p{A)\ 



n— too q 
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Once this is established, since sets of the form [—00, a) are convergence determining, we will 
have weak convergence. We proceed by dividing into three terms using the triangle inequal- 
ity and the set-algebra identity that \fnA + pn = (y/rA + pr) (J {y/ra + pr, yfna + pn]: 



'T n - pn 



n 



€A\- p(A) 



(/%) (y/nA + pn) - p{A)\ 
Co.rA + ^2 C r)n jl r (y/nA + pn) - p(A) 



< 



r=2 
n-Bj 



+ 



+ 



C ,n<5o + ^2 C^nfir [y/nA + pn 
r=2 J 
\ 

C r,nPr [y/rA + pr) - p(A) 



y r=n-d n +l 
n 



C r>n p, r \ [y/ra + pr, y/ria+pii\ 

r=n-d n +l 



We now show that each term individual goes to zero as n — > 00. Each term is handled in a 
separate claim. 

Claim. 



Co >n So + ^ C r ,nPr (VnA + pn) 



r=2 







Proof of Claim. This is immediate from the second property in 14.41 using the following 
relation. 



Co,n^O + ^ C r,nPr [y/nA + 



pn) 



r=2 



< 



r=2 



Claim. 



□ 



C r ^ n p T [yftA + pr) - p{A) 
K r=n-e n + l J 

Proof of Claim. This follows from the results of 14.41 and from the weak convergence for T n , 
namely Tn /^ n P- Fix e > 0. We use the weak convergence of T n applied to the set A to 
find an N G N so large so that \p r (y/rA + pr) - p(A)\ < e for all r > N — 9n + 1. Then, 
for all n > N we have that: 
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C r,nVr {VrA + pr) - p(A) 

r=n—9 n +l 



< 



°r,n (fir (V^A + pr) - p(A)) 

r=n—9„+l 



+ 



< 



I n—On \ 

\C ,n + J2 Cr A P^ A ) 

n \ / n-6>„ \ 

^] C r ,n J e + I Co, n + Yj Cr,n J 
=n-0 n +l / V r=2 / 



1 • e + • p(^) = e 



where the limits are from 14.41 Since this tends to e as n — > oo, it is eventually less than 2e 
for sufficiently large n. Since this holds for all e > 0, we have the result of the claim. □ 

Claim. 



C r}Tl fl r (y/ra + pr, y/na + pv/ | 

. r=n-0 n +l 







Proof of Claim. Fix an e > 0. Since the measure p has no atoms, choose a 5 > so small so 
that p(a,a + 5] < e. Now, since Tn //^ n p, choose N so large so that for all n > N — 9jq + 1 
we have: 



T n -pn 



n 



G (a, a + 5] — p(a,a + 5] 



< € 



Now choose M so large, so that the following holds whenever n > M and n > r > n — 6 n + l: 

\fn — \fr n — r 
a- 



y/r v r 

Such an M always exists because the inequalities n > r > n — 8 n + 1 give: 



>n — y/r n — r 
a- ^ hp — r=- 



n \ n — r 
r / Jr 



< a| 



1 



1 +P 



n 



1 /X 



'2 V n 



So we see that the limit works because of our choice of n so that 9 n satisfies ^= — > 0. 

For n larger than both iV and M then, we will have then for every r with n > r > 
n — 9 n + 1 that: 



p r (y/r a + pr, yna + pre] 



< P 



T r — pr 



T r — pr 



G a,a+ a 



n — \ r n — r 



+ p- 



G (a, a + 5] 



< p (a, a + 5] + e 

< e + e = 2e 
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Hence, for such n we have: | (Ylr=n~e n +i c r,np-r) (\/™ + pr, y/na + pn\\ < 1 ■ 2e. Since this 
holds for every e > 0, we get the desired result. 
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